A tight-binding approach to uniaxial strain in graphene 
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We analyze the effect of tensional strain in the electronic structure of graphene. In the absence 
of electron-electron interactions, within linear elasticity theory, and a tight-binding approach, we 
observe that strain can generate a bulk spectral gap. However this gap is critical, requiring threshold 
deformations in excess of 20%, and only along preferred directions with respect to the underlying 
lattice. The gap less Dirac spectrum is robust for small and moderate deformations, and the gap 
appears as a consequence of the merging of the two inequivalent Dirac points, only under considerable 
deformations of the lattice. We discuss how strain-induced anisotropy and local deformations can 
be used as a means to affect transport characteristics and pinch off current flow in graphene devices. 
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It is now well established that sp^ bonded carbon sys- 
tems feature record-breaking mechanical strength and 
stiffness. Investigations in the context of carbon nan- 
otubes reveal intrinsic strengths^ that make these sys- 
tems the strongest in nature, Recently, graphene — the 
mother of all sp'^ carbon structures — has been confirmed 
as the strongest material ever to be measured^, being 
able to sustain reversible elastic deformations in excess 
of 20%\ 

These mechanical measurements arise at a time where 
graphene draws considerable attention on account of 
its unusual and rich electronic properties. Besides the 
great crystalline quality, high mobility and resilience 
to high current densities^, they include a strong field 
effect^, absence of backscattering^ and a minimum metal- 
lic conductivity''. While many such properties might 
prove instrumental if graphene is to be used in future 
technological applications in the ever pressing demand 
for miniaturization in electronics, the latter is actually 
a strong deterrent: it hinders the pinching off of the 
charge flow and the creation of quantum point contacts. 
In addition, graphene has a gapless spectrum with lin- 
early dispersing, Dirac-like, excitations^'^. Although a 
gap can be induced by means of quantum confinement 
in the form of nanoribbons^^ and quantum dots^^, these 
"paper-cutting" techniques are prone to edge roughness, 
which has detrimental effects on the electronic properties. 
Hence, a route to induce a robust, clean^ bulk spectral 
gap in graphene is still much in wanting. 

In this paper we inquire whether the seemingly inde- 
pendent aspects of mechanical response and electronic 
properties can be brought together with profit in the con- 
text of a tunable electronic structure. Motivated by re- 
cent experiments showing that reversible and controlled 
strain can be produced in graphene with measurable 
effects^^'^^, we theoretically explore the effect of strain 
in the electronic structure of graphene within a tight- 
binding approach. Our calculations show that, in the 



absence of electron-electron interactions, a gap can be 
opened in a pure tight binding model of graphene for 
deformations beyond 20%. This gap opening is not a 
consequence of a broken sublattice symmetry but due to 
level crossing. The magnitude of this effect depends on 
the direction of applied tension, so that strain along a 
zig-zag direction is most effective in overcoming the gap 
threshold, whereas deformations along an armchair di- 
rection do not induce a gap. Unfortunately, such large 
threshold deformations render strain an ineffective means 
to achieve a bulk gapped spectrum in graphene. We dis- 
cuss alternate means to impact transport and electronic 
structure using local strain profiles. 



I. MODEL 

We consider that electron dynamics of electrons hop- 
ping in the honeycomb lattice is governed by the nearest 
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FIG. 1: (Color online) (a) Tension geometry considered in 
the text. The zig-zag direction of the honeycomb lattice is 
always parallel to the axis Ox. 
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the stress and strain tensors. In general the components 
Cijki depend on the particular choice of the Cartesian 
axes. Incidentally, for an hexagonal system under pla- 
nar stress in the basal plane, the elastic components are 
independent of the coordinate system. This means that 
graphene is elastically isotropic^^. 

The analysis of strain is straightforward in the princi- 
pal system Ox'y' where we simply have T = Tcx' ' 
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^ijkl^kl ~ T Sijkl^kx^lx — T S', 
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a) Honeycomb lattice geometry. The 
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FIG. 2: (Color online) 

vectors Si — a(^,— |), 82 — a(0, 1) 5^ — a(— |j con- 
nect A-sites (red/dark) to their B-site (blue/light) neighbors, 
(b) The first Brillouin zone of undeformed graphene, with its 
points of high symmetry. 



Given that only five compliances are independent in 

graphite (viz., Sxxyyi ^xxyyi ^xxzzi ^zzzzi ^yzyz^ 1 1^ 

follows that the only non-zero deformations are 
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neighbor tight-binding Hamiltonian 

't{R,S)a'^{R)b{R^S)- 



which represent the longitudinal deformation and Pois- 
son's transverse contraction. If we designate the tensile 



strain hy e = T Sx 



the strain tensor can be written 



H 



H. c. . (1) terms of Poisson's ratio, a = —Sxxxy/Sx 



Here R denotes a position on the Bravais lattice, and 5 
connects the site R to its neighbors; a{R) and h{R) are 
the field operators in sublattices A and B. The first thing 
to emphasize is that, under general stress conditions, the 
hopping t{R^ 5) will be generally different among differ- 
ent neighbors. We are interested in the elastic response, 
for which deformations are affine. This means that even 
though the hoppings from a given atom to its neighbors 
can be all different, they will be the same for every atom. 
Therefore, as depicted in Fig. 2(a), we need only to con- 
sider three distinct hoppings: ti = t(5i), ^2 = ^(^2), and 
ts = t{5^). The relaxed equilibrium value for ti = t{6i) 
is to ^ 2.7 eV^. Our goal is to investigate the changes 
that strain induces in these hoppings, and what impact 
they have in the resulting electronic structure. 

Throughout this paper we shall use the C-C equilib- 
rium distance, ao = 1.42 A, as unit of length, and will 
frequently use to as unit of energy. 
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This form shows that graphene responds as an isotropic 
elastic medium. For Poisson's ratio we use the value 
known for graphite: a = 0.165^^. It should be mentioned 
that when stress is induced in graphene by mechanically 
acting on the substrate (i.e. when graphene is adhering 
to the top of a substrate and the latter is put under ten- 
sion, as is done in Ref. 12), the relevant parameter is in 
fact the tensile strain, rather than the tension T^^. For 
this reason, we treat e as the tunable parameter. Since 
the lattice is oriented with respect to the axes Oxy, the 
stress tensor needs to be rotated to extract information 
about bond deformations. The strain tensor in the lattice 
coordinate system reads 



cos — a sin 6 (1 + <j) cos^sin6> 
{1 -\- a) cos 6 sin 6 sin^ 6 — a cos'^ 6 



(6) 



II. ANALYSIS OF STRAIN 

We are interested in uniform planar tension situations, 
like the one illustrated in Fig. 1(a): the graphene sheet is 
uniformly stretched (or compressed) along a prescribed 
direction. The fixed Cartesian system is chosen in a way 
that Ox always coincides with the zig-zag direction of 
the lattice. In these coordinates the tension, T, reads 
T = T cos{0) Cx^T sin{e) Cy . 

As for any solid, the generalized Hooke's law relating 
stress, Tij and strain Sij has the form 

^ij — Cijkl Ski •) Sij = Sijki Tj^i (2) 

where Cijki {Sijki) are the components of the stiffness 
(compliance) tensor. Since we address only states of pla- 
nar stress, we resort to the 2-dimensional reduction of 



III. BOND DEFORMATIONS 

If represents a general vector in the undeformed 
graphene plane, its deformed counterpart is given, to 
leading order, by the transformation 

v = {l^s)-v^. (7) 

Especially important are the deformations of the nearest- 
neighbor bond distances. Knowing Sij one readily ob- 
tains the deformed bond vectors using (7). The new bond 
lengths are then given by 

l^il = 1 + 1^11-^^12 + 1^22 (8a) 
1^2! = 1 + 6:22 (8b) 

l^sl = 1 + 1^11 + ^^12 + 1^22 (8c) 
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FIG. 3: (Color online) (a) Plot of ti/t2 vs t3/t2 as a func- 
tion of strain, and 0. Closed lines are iso-strain curves, 
and arrowed lines correspond to the trajectory of the point 
(ti/t2, t3/t2) as s increases, calculated at constant angle. The 
graph is symmetric under reflection on both axes. In the 
shaded area the spectrum is gapless. The blue iso-strain line 
{e ^ 0.23) corresponds to the gap threshold. In panel (b) 
we show the angular dependence of ti,2,3 for e = 0.05 and 
£ = 0.23. 



Of particular interest are the cases ^ = and 6 = 7r/2 
since they correspond to tension along the zig-zag (Z) 
and armchair {A) directions: 



Z : 
A: 



l^l| = l^3| = l' 
l^l| = l^3| = l' 



ea . 



ea . 
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The modification of these distances distorts the re- 
ciprocal lattice as well, and the positions of the high- 
symmetry points shown in Fig. 2(b) are shifted. The 
primitive vectors of the reciprocal space are denoted by 
6i^2, and, in leading order, change according to 
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Most importantly, the symmetry point K = (^^,0) 

(that coincides with the Fermi point in the undoped, 
equilibrium situation, and chosen here for definiteness) 
moves to the new position 
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for a general deformation, and in leading order in strain. 
For uniaxial tension this reduces to 



K 



47r 
3V^ 



1 



£(l-a). 



-6:(1 + cr) sin[2l9] 



(12) 



The factor of 20 means that the shift is the same for 
the A and Z directions, in leading order. These general 
results will be important for our subsequent discussion. 



IV. HOPPING RENORMALIZATION 

The change in bond lengths (8) leads to different hop- 
ping amplitudes among neighboring sites. In the Slater- 
Koster scheme^ ^, the new hoppings can be obtained from 
the dependence of the integral Vppj^ on the inter-orbital 
distance. Unfortunately determining such dependence 
with accuracy is not a trivial matter. Many authors re- 
sort to Harrison's flyleaf expression which suggests that 
VppT^{l) oc However this is questionable, inso- 

far as such dependence is meaningful only in matching 
the tight-binding and free electron dispersions of simple 
systems in equilibrium (beyond the equilibrium distance 
such dependence is unwarranted^^). It is indeed known 
that such functional form fails away from the equilibrium 
distance^^, and a more reasonable assumption is an ex- 
ponential decay In line with this we assume that, in 
graphene 



VppTr{l) = toe 



-3.37{l/ao-l) 



(13) 



where the rate of decay is extracted from the experi- 
mental result dVppT^/dl = —6.4 eV / A?^ . As a consis- 
tency check, we point out that, according to Eq. (13), 
the next-nearest neighbor hopping {t') would have the 
value VppT^{\/?>aQ) = 0.23 eV, which tallies with existing 
estimates of f in graphene^. 



V. GAP THRESHOLD 

The bandstructure of Eq. (1) with arbitrary hoppings 
^1,^2, ts is given by 



E{kx,ky) = ± \t2 + fee' 
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(14) 



Here both and the primitive vectors (see Fig. 2(a) 
for the definition of the vectors a^) change under strain: 
the hoppings change as per Eqs. (13,8), and the primitive 
vectors as per Eq. (7). This generalized dispersion has 
been previously discussed in Refs. 22,23, under the as- 
sumption that only the hopping elements change, with- 
out lattice deformation. It was found that the gapless 
spectrum is robust, and that a gap can only appear un- 
der anisotropy in excess of 100% in one of the hoppings. 
More specifically, the spectrum remains gapless as long 
as the condition 
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(15) 



is in effect. This condition corresponds to the shaded 
area in Fig. 3(a). Using the results in eqs. (6, 8, 13) we 
have mapped the evolution of the hoppings with e and 0. 
This allows us to identify the range of parameters that 
violate (15), and to obtain the threshold for gap open- 
ing. For a given 0, we follow the trajectory of the point 
(ti/fe, fe/fe) as strain grows, starting from the isotropic 
point at £ = 0. The result is one of the arrowed curves in 
Fig. 3(a). The value of e at which this curve leaves the 
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FIG. 4: (Color online) Top row shows density plots of the 
energy dispersion, E{kx,ky), for {e = 0, ^ = 0} (a), {e — 
0.2, = 7r/2} (b), and {e = 0.2, 6* = 0} (c). In these plots, 
the central white dashed lines represent the boundary of the 
first BZ of the undeformed lattice, while the solid white lines 
mark the boundaries of the BZ for the deformed lattice. In (d) 
we have a cut of (c) along /c^ = 0, showing the merging of the 
Dirac cones as strain increases, and the ultimate appearance 
of the gap. In panel (e) we compare the gap given by Eq. (18) 
(line) with the result obtained from direct minimization of the 
energy in the full BZ (dots). 



shaded area corresponds to the gap threshold for that 
particular angle. From such procedure, summarized in 
Fig. 3(a), we conclude that: (i) the gap threshold is at 
e ^ 0.23 (~ 20%); (ii) the behavior of the system is 
periodic in 6 with period 7r/3, in accord with the symme- 
try of the lattice; (iii) tension along the zig-zag direction 
{0 = 0, 7r/3, . . . ) is more effective in overcoming the gap 
threshold; (iv) tension along the armchair direction never 
generates a gap. 

The two panels of Fig. 3(b) contain plots of the indi- 
vidual ta for two particular values of strain. It is clear 
that, for deformations along the Z direction, the highest 
relative change occurs along the zig-zag bonds (^1,3), and 
conversely for deformations along the A direction. This 
could have also been anticipated from Eqs. (9) and the 
smallness of a. 



VI. CRITICAL GAP 

The fact that the isotropic point (1, 1) in Fig. 3(a) is 
surrounded by an appreciable shaded area, means that 
the gapless situation is robust, and the emergence of the 
gap requires a critical strain. The physical effect be- 
hind such critical gap lies in the fact that, under strain, 
the Dirac cones drift away from the points K' in the 



Brillouin zone (BZ). Before we proceed, it is pertinent to 
advance a crucial detail: in a deformed lattice, the Dirac 
points (i.e. the positions in the BZ where conduction and 
valence bands touch conically) and the symmetry points 
Ki do not coincide. In what follows we shall distinguish 
them explicitly. 

To be more definite, we examine the position of the 
minimum energy^^ for the bands obtained from Eq. (14), 
which can be done exactly if we assume that the lattice 
remains undeformed. Due to the particle-hole symmetry, 
we minimize E{k)'^. Let us assume that ti = ts 7^ ^2, 
which applies for tension along either zig-zag or arm- 
chair directions. In that case the spectrum has minima 
at exactly 
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and all symmetry related points. The ± sign refers to 
one possible choice for the two inequivalent valleys. In 
addition to these local minima, that correspond to the 
Dirac points, the dispersion has saddle points at 
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(and all symmetry related points). These are just the 
points TVfi, 71^2, Ms shown in Fig. 2(b) and their posi- 
tion is independent of U'^^. The values of energy at these 
points aTeE{Mi) = l^i+ts-tal, E{M2) = l-^i+ts+tal 
and E{Ms) = \ti-t2^t3\. 

The result in Eq. (16) shows that the Dirac points drift 
away from the K point, and the direction of that drift 
is dictated by the relative variations in ti. For example, 
for uniaxial tension along the Z we have ^2 > ti = ^3, 
and therefore the minimum of energy moves to the right 
(left) of Ki {K{) [cfr. Fig. 2(b)]. This means that the 
inequivalent Dirac points move toward each other, and 
will clearly meet when 2ti = ^2- They meet precisely at 
the position of the saddle point M2. Throughout this 
process, the dispersion remains linear along the two or- 
thogonal directions, albeit with different Fermi veloci- 
ties. If the hoppings change further so that 2ti > ^2, 
the solution (16) is no longer valid, and the minimum 
lies always at M2. Since the energy at this saddle point 
is given exactly by E{M2) = \2ti — ^2!, the system be- 
comes gapped, with a gap A = 2|2ti — t2\. Moreover, 
the dispersion becomes peculiar in that it remains linear 
along one direction (the y direction in this example) and 
quadratic along the other. The topological structure is 
also modified since the two inequivalent Dirac cones have 
merged^^. 

These considerations assume that the hoppings can 
change but the lattice remains undeformed. Under a real 
deformation both lattice and hoppings are affected. The 
lattice deformation will distort the BZ but will not af- 
fect the aspects discussed above. This can be clearly 
seen from inspection of the energy dispersions plotted in 
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FIG. 5: (Color online) Closeup of the energy dispersion E(k) 
in the vicinity of the Dirac points. Panel (a) shows a pre- 
critical situation (in which the two cones persist) for a strain 
of s = 0.2 along the zig-zag {0 = 0) direction. Notice how the 
saddle point approaches zero energy with increasing strain. 
In (b) we show the bandstructure at precisely the critical 
strain for which the two Dirac cones meet. The dispersion is 
quadratic along kx (the direction of strain) but remains linear 
along ky. This can be seen in (c) where we show the critical 
dispersion viewed from the kx and ky axes, respectively. 



Fig. 4(a-c). The plotted dispersions include the deforma- 
tion of the BZ and the change in the hoppings, simulta- 
neously. For strain along the A direction the nonequiv- 
alent Dirac cones move in opposite directions and never 
meet [Fig. 4(b)]. However, if the deformation is along 
the Z direction, the cones always approach each other 
[Fig. 4(c)], and will eventually merge. This merging is 
seen in detail in Fig. 4(d) where a cut along ky = is 
presented to show the emergence of the gap beyond the 
threshold deformation. For tension along an arbitrary 
direction (except armchair) the cones always merge, the 
>Zdirection being the optimal orientation, requiring less 
strain [cfr. Fig. 3(a)]. Precisely at the critical point, the 
dispersion is linear along kx and quadratic along ky^ as 
shown in Fig. 5. This modification in the dispersion along 
one of the Cartesian directions has peculiar implications 
for the DOS and Landau level quantization^^. 

The gap is a result of this Dirac cone merging process, 
and the origin of the high critical strain is now clear: one 
needs to deform enough to bring the two Dirac points to 
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FIG. 6: (Color online) (a) A close-up of the energy dispersion 
close to its minimum, for tension along ^ = and e — 0.05. 
The solid white lines show the intersection of the Bragg planes 
that define the boundary of the first BZ in the deformed lat- 
tice, while the dashed white lines represent the same boundary 
es in equilibrium. It is clear that the Dirac point lies neither 
at K ^ nor at its deformed counterpart. The energy contours 
are labeled in eV. (b) A contour plot of the absolute value of 
the Fermi velocity, hvp — VkE{k) for the same region shown 
in (a). 

coincidence. This agrees with the existing understanding 
that the gapless Dirac spectrum in graphene is robust 
(topologically protected) with respect to small perturba- 
tions. 

For strain along ^ = 0, as discussed above, the gap is 
conveniently given by 

Eg{e) = 2 |2ti(£) - t2{e)\e{t2 - 2ti) . (18) 

An example of the strain dependence of Eg can be seen 
in Fig. 4(e). In it we see the agreement between the gap 
given by Eq. (18) and the value extracted from a direct 
minimization of E(kx-,ky) in the full (deformed) BZ. 

From Fig. 4(b) one can see that pulling along an arm- 
chair direction imparts ID-like features to the system: 
the dispersion becomes highly anisotropic. This is ex- 
plained on account of the results plotted in Fig. 3(b) 
which show that stress along A tends to weaken one 
bond only. In extreme cases, the weak bond can be 
highly suppressed leaving only a set of ID chains^^. This 
means that strain along certain directions can be used as 
a means to induce preferred anisotropy in electric trans- 
port. In contrast, pulling along a zig-zag direction tends 
to dimerize the system for large deformations, which ul- 
timately explains the appearance of the gap in this case. 

VII. POSITION OF THE DIRAC POINT 

The fact that there are two concurrent effects deter- 
mining the changes in the bandstructure (viz. the lat- 
tice distortion itself, and the modification in the nearest- 
neighbor hoppings) means that the position of the min- 
imum of energy does not coincide with the symme- 
try points of the deformed BZ. This is documented in 
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Fig. 6(a) where we provide a close-up of the energy dis- 
persion close to the Dirac point. This should be clear 
from the foregoing discussion on the merging of the Dirac 
points. In any case we want to stress this effect and illus- 
trate it by analyzing small perturbations with respect to 
the undeformed situation. For definiteness let us focus 
again in the case ti = t^. The position of the new K was 
given already in eqs. (11,12). The position of the Dirac 
point is given by eq. (16) when only the hoppings change, 
but not the lattice. For a small perturbation (^2 ~ ti) 
the result (16) reduces to^^ 
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(19) 



One can calculate the correction to this result simulta- 
neously accounting for the lattice deformation. But the 
lengthy expression that results is less important than the 
qualitative effect: the corrections to the expression (19) 
depend on the specific details of the variation of ti with 
distance. Consequently, the Dirac point and the K point 
of the deformed lattice do not coincide in general. The 
equilibrium situation, in which they coincide, is a very 
particular case. 

In fact, even assuming a simple lattice distortion that 
does not change the hoppings will move the Dirac point 
away from the symmetry point of the resulting lattice. 
This can be seen from a low energy expansion of (14) 
putting ti = t. The position of the Dirac point Ki in 
such case moves to 
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which is clearly different from (11). This fact is of critical 
relevance when interpreting results of similar calculations 
obtained ab-initio, as will be discussed below. 

The fact that the Dirac point drifts from the corner 
of the BZ means that there are no longer 3 equivalent 
pairs of points in the first BZ for the neutral system, 
but only one pair of non-equivalent points in general (in 
other words, in the undoped and undeformed lattice the 
Fermi surface is distributed among the 6 degenerate K 
points in the boundary of the BZ, whereas for a general 
deformation we have only two within the rst BZ). In the 
situation shown in Fig. 6(a), for example, the Dirac point 
shown in the figure lies outside the first Brillouin zone. 
The one inside the first BZ is actually the (equivalent) 
Dirac point that moved away from K2 (or i^s), in the 
notation of Fig. 2(b). 



VIII. DISCUSSION 

We have seen that, within the tight-binding Hamilto- 
nian written in Eq. (1), uniform tension can induce a 
bulk spectral gap in graphene. However, at least within 
a non-interacting tight-binding approach, the gap thresh- 
old is very difficult to overcome, if at all possible. Since 



a tensional strain in excess of 20% is required to observe 
such feature, several comments are in order. 



A. On the approximations employed 

We start by noticing that in our calculation we kept 
only the lowest order terms in s. In addition, al- 
though strain magnitudes of ~ 20% are not unreason- 
able, graphene is expected to be in the non-linear elastic 
regime at those deformations^. Therefore, non-linear cor- 
rections can be relevant at the quantitative level in the 
vicinity of the threshold. 

Notwithstanding, our main result is robust: no gap 
can be opened under planar tension situations, except 
in highly strained situations. This conclusion does not 
depend on having taken a linear approximation insofar 
as it should be valid up to deformations in the range of 
5-10%. 

With respect to our tight-binding parametrization in- 
cluding only nearest neighbor hopping, we should men- 
tion that Kishigi et a/. have shown that inclusion of next- 
nearest neighbor terms {t') can, alone, generate a gap^^. 
But this requires a very specific deformation of the lat- 
tice, unlikely to occur under simple tension. The presence 
of f can also lead to other effects, like tilted Dirac cones 
as discussed in Ref. 31. 

It is expected that the planar arrangement of carbon 
atoms in freely hanging graphene should become unstable 
with respect to a buckled or rippled configuration, or 
even experience mechanical failure for moderate to high 
tension. The presence of a substrate should provide more 
stability for the planar distribution of the carbon atoms. 
In fact, a recent experiment published during revision 
of this manuscript^^ shows that reversible strain of the 
order of 18% can be induced in graphene deposited on 
flexible plastic substrates. 



B. On related ab-initio calculations 

Secondly, some ab-initio calculations seem to show 
that a gap is present in graphene for arbitrarily small 
tensions^^'^^. But these reports have some conflicting 
details. For example there is an order of magnitude dis- 
crepancy between the gap predicted in these two refer- 
ences for 1% strain. In addition, Ref. 33 claims their 
ab-initio result agrees with the bandstructure (14) after 
a suitable choice of hoppings. As we showed above this 
cannot be the case, since there is always a (large) thresh- 
old for the appearance of the gap. Consequently, further 
clarification regarding ab-initio under strain is desired. 

One issue that requires special attention in inter- 
preting density functional theory (DFT) calculations of 
graphene's bandstructure under strain is the shift of the 
Dirac point. As we stressed earlier, when graphene is 
strained the Dirac point (position of the energy mini- 
mum) does not lie at any symmetry point of the lat- 
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tice. This is paramount because DFT calculations of 
the bandstructure rely on a preexisting mesh in recipro- 
cal space, at whose points the bandstructure is sampled. 
These meshes normally include points along high sym- 
metry lines of the BZ. But in the current problem, the 
use of such a traditional mesh will not be particularly 
useful to distinguish between a gapped and gappless sit- 
uation. Since the Fermi surface of undoped graphene is a 
single point, unless the sampling mesh includes that pre- 
cise point, one will always obtain a gap in the resulting 
bandstructure. 

Most recently, we became aware of two new devel- 
opments from the ab-initio front, that shed the light 
needed to interpret the earlier calculations mentioned 
above. One of them^^ is a revision of the DFT calcu- 
lations presented in Ref. 12. In this latter work it is 
shown that, upon careful analysis, the DFT calculation 
shows no gap under uniaxial deformations up to ~ 20%. 
In fact the authors mention explicitly that the shift of 
the Dirac point from the high symmetry point mislead 
the authors in their initial interpretation of the band- 
structure. In another preprint^^, independent authors 
show that their DFT calculations reveal, again, no gap 
for deformations of the order of 10% in either the ^ or ^ 
directions. These subsequent developments confirm our 
prediction that only excessive planar deformations are 
able to engender a bulk spectral gap in graphene. 



C. Anisotropic Transport 

Several effects of tensional strain are clear from our re- 
sults. Tension leads to one-dimensionalization of trans- 
port in graphene by weakening preferential bonds: trans- 
port should certainly be anisotropic, even for small ten- 
sions. One example of that is seen in Fig. 6(a) where a 
strain of 5% visibly deforms the Fermi surface. Fig. 6(b), 
where the Fermi velocity is plotted for the same region in 
the BZ, further shows that, for the chosen tension direc- 
tion, the Fermi surface is not quite elliptical but slightly 
oval, in a reminiscence of the trigonal warping effects. 

The anisotropy in the Fermi velocity can become quite 
large, as shown in Fig. 7. There we plot the ratio 
'^F(max)/'^F(min) foi" = 50 mcV. As cau be sccu in 
the cuts of Fig. 4(d), near the critical strain the Fermi 
level may touch the van Hove singularity (midpoint be- 
tween the two cones) where vp = along one direction, 
which leads to formal divergence of the ratio between the 
maximum and minimum Fermi velocities. 

The Fermi surface anisotropy has been captured 
in transport experiments that reveal a considerable 
anisotropy in the resistivities parallel (Rxx) and per- 
pendicular (Ryy) to the tension direction^^. The au- 
thors of Reference 32 report in their Fig. 4c resistivity 
anisotropics of up to one order of magnitude at ~ 19% 
strain. We can make a simple estimate of the anisotropics 
expected in the light of our results for the deformed band- 
structure. For that we follow a Boltzmann approach to 



the relaxation time and DC conductivity. As is usually 
the case, the relevant electron states are the ones lying 
in a narrow vicinity of the Fermi surface. On account 
of the Fermi surface anisotropy, we can take the Fermi 
velocities in the direction of the electric field only. For 
example, in the situation shown in Fig. 6(a) for uniaxial 
strain along Ox, the relevant vp^s determining R^x and 
Ryy will be the ones along the major and minor axes, 
respectively. 

For the purposes of our estimate, we take the Boltz- 
mann longitudinal conductivity for scattering out of un- 
screened charged impurities of valence Z and concentra- 
tion ni^^, given by 



a = 



2e^ Trh'^VpTi 
h Z'^e^rij 



(21) 



and follow the reasoning above: replacing vp Vp-Up^ 
where Up is the direction of the electric field in each case 
(parallel or perpendicular to the tension direction). We 
immediately see that this estimate leads to a resistance 
anisotropy of^*" 



Rxx 

R~ 



^F(max) 
yy \ ^F(min) 



(22) 



For the maximum strain used in the experiment (~ 19%) 
our plot in Fig. 6(b) yields a ratio of 2.8, which, per (22), 
corresponds to a resistance anisotropy of roughly 8 fold, 
consistent with the measured anisotropy. 

Our results apply to exfoliated and epitaxial graphene 
alike. As matters currently stand, it is perhaps more rel- 
evant in the context of the latter, since graphene grown 
epitaxialy on SiC is almost always under strain^^. The 
strained configurations are imposed by the lattice mis- 
match with the substrate, and can be controlled by 
changing the growth and annealing conditions'^. For 
these systems, the relaxed starting configuration is al- 
ready deformed. 

Lastly, it is important to point out that, even though 
a spectral gap seems to require extreme strain, one can 
generate a transport gap by means of local, small, defor- 
mations. It has been shown in Refs. 40,41 that tunneling 
across a locally strained region is highly suppressed, and 
leads to a transport gap (i.e., a suppression of electrical 
conductivity) at small densities, even in the absence of a 
bulk spectral gap. 



IX. CONCLUSIONS 

Within a non-interacting, nearest-neighbor tight- 
binding approach, we have shown that opening a spec- 
tral gap in strained graphene requires deformations of 
the order of 20%. This result is confirmed by the most 
recent ab-initio calculations. Such an extreme strain is 
required on account of the stability of the Dirac points 
in graphene, that renders the spectrum gappless unless 
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FIG. 7: (Color online) Anisotropy in the Fermi velocity, as a 
function of strain for a Fermi energy of 50meV. The vertical 
axis shows the ratio of the maximum to the minimum values 
oi vf in the entire Fermi surface. 

the two inequivalent Dirac points merge. The merg- 
ing requires substantial anisotropy in the hopping inte- 



grals, only achieved under high strain. General features 
of strained graphene are an anisotropic Fermi surface, 
anisotropic Fermi velocities, and the drift of the Dirac 
points away from the high symmetry points of the lat- 
tice. 

Uniform planar stain appears to be an unlikely can- 
didate to induce a bulk gap in graphene. Nevertheless, 
strain (local or uniform) can be an effective means of tun- 
ing the electronic structure and transport characteristics 
of graphene devices. Even if the bulk gap turns out to 
be challenging in practice, local strain could be used as 
a way to mechanically pinch off current flow. 
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